Phase- dependent noise correlations in normal-superconducting structures 
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We study nonequilibrium noise correlations in diffusive normal-superconducting structures in the 
presence of a supercurrent. We present a parametrization for the quasiclassical Green's function in 
the first order of the counting field x- This we employ to obtain the voltage and phase dependence 
of cross and autocorrelations and to describe the role played by the setup geometry. We find that 
the low- voltage behavior of the effective charge ^eff describing shot noise is a result of a competition 
between anticorrelation of Andreev pairs due to proximity effect and the depression of the local 
density of states. Furthermore, we show that the noise correlations are independent of the sign of 
the supercurrent. 6 
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The charge transmitted through a disordered conduc- 
tor in unit time varies due to the quantum nature of 
the transport process. This variation is characterized by 
the current distribution, whose width at low tempera- 
tures is directly related to the shot noise. ^ In metallic 
conductors, the corrections induced by the quantum co- 
herence on the current and conductance distributions are 
small^i^ but in a metal in contact to a superconductor, 
more pronounced effects may be observed, e.g., in the 
out-of-equilibrium noise experiments j^^^i^i^ In the context 
of a normal-superconducting (NS) two-terminal setup, 
the theoretical endeavours have recently covered, e.g., the 
voltage dependence of the shot noise^^^ and certain low- 
bias anomalies. ^^^^ Multiterminal structures have been 
discussed in the incoherent regime, ^^^^^^^^^^^ and in the 
presence of a supercurrent, in a short junction^ and for 
specific values of a phase difference in a three-terminal 
setup. The latter was described by a method based on 
a direct discretization of the equations governing the full 
counting statistics [see Eq. ([I])].^ 

In the presence of the superconducting proximity effect 
and at voltages of the order of Et/c^ the effective charge 
geff characterizing the magnitude of shot noise is lower 
than 2e, the value corresponding to a Cooper pairi^i^i^ 
Here Et = hD / is the Thouless energy, with D the 
diffusion constant and L the wire length. Applying a 
supercurrent in a three-terminal structure (Fig. [1]) allows 
one to study the nature of this proximity-induced change 
in g'eff- For example, with this approach, we directly 
show that the lowering of q^^ is due to a competition 
of anticorrelation effects induced by the superconducting 
proximity effect^ and the depression of the local density 
of states. In principle, the effect of an extra current 
in the system would be twofold: to tune the coherent 
effects and to induce its own correlations. Here we show 
that within our noninteracting model supercurrent does 
only the previous as all the correlations in the system are 
independent of the supercurrent sign. 

The full counting statistics^-ii^i^ has recently become 
the method of choice to calculate shot noise in diffu- 
sive mesoscopic conductors but, to our knowledge, the 
cross correlations in the presence of supercurrent have 




FIG. 1: Setup schematic studied in this work. 



not been studied. The statistics is accessed by perform- 
ing a counting rotation of the Green's function qq in one 
of the terminals, using the counting field x-^^ To access 
the noise correlations, one has to expand the resulting 
function g{x) = 9o — ^(x/2)^i + O(x^) in powers of the 
counting field and solve the problem in the first order in 
X' Within quasiclassical formalism, the spectral quanti- 
ties related to go may be represented by two parameters, 
and (j) characterizing the magnitude and phase of the 
pair correlations. The distributions of the electrons and 
holes are treated by dividing the functions into even and 
odd components with respect to the Fermi surface, fr 
and /l-^^ In the absence of supercurrent, a parametriza- 
tion for gi has been given in Ref. [TO. Here we present a 
parametrization for gi applicable also for a finite super- 
current. Since the resulting essentially linear method is 
considerably faster than the previous "direct discretiza- 
tion" approach, we have been able to extensively study 
the roles played by the setup geometry, dissipative cur- 
rents, coherence and phase gradients. 

In our setup (Fig. [1]), the two superconducting termi- 
nals Ti and T2 are connected by diffusive wires 1 and 2 
and, at the central node, by a diffusive control wire 3 to a 
normal reservoir, T3. The phase difference Acj) = (1)2 — 4>i 
between the superconductors may be generated by fabri- 
cating a superconducting loop and applying a magnetic 
flux through the loop or by an external driving of super- 
current. The lengths of the wires are denoted by ^^1,2,3 
and the currents into the terminals /i,2,3- We designate 
the cross section of the control wire by Ac and suppose 
that the other wires are equally wide, with cross section 
A. The electric potential is assumed to vanish in the 
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superconductors. We assume good contacts at the inter- 
faces and a vanishing temperature and consider voltages 
below the superconducting energy gap V <C A/e. 

In the quasiclassical diffusive limit, the triangular 
matrix go = g{x = 0) in Nambu(")-Keldysh(~) space 
may be expressed through A, and which may 
be parametrized using A = —t^W%^ K = Rh — hA, 
h = fL ^ /tT3, R = cosh(6>)f3 H- sinh(6>)(cos(^)if2 + 
sin(0)ifi).^^ Here 6 characterizes the strength of the su- 
perconducting proximity effect, (j) is the superconducting 
phase, and are the Pauli matrices in Nambu space. The 
counting field which we introduce in the normal reser- 
voir, appears as a gauge transformation of the Green's 
function in the same terminal. At a vanishing tempera- 
ture, it suffices to concentrate only on the energy regime 
< e < eV. In this case, without counting rotation, 
one has in the normal terminal gN,o = 9n{x = 0) = 
T3 0-3 + To (o^i + ^0-2), where are Pauli matrices in 
Keldysh space. In the superconducting terminals we have 
gs = (cos(0)f2 + sin((/))fi) ^ao. The generalized Green's 
functions g{x) obey the Usadel equation^^ similar to that 
in the zeroth order in x 

- J^^^^^^) = -ie[f^,g{x)] , (1) 

with J{x) = —LGDg{x)dxg{x)^ fs = % ^a^. Here Gd is 
the normal-state conductance of the wire with length L, 
X is the coordinate along the wire, and e is the energy. 
The Green's function satisfies the normalization condi- 
tion g'^{x) = To ao. At the NS interfaces the Nazarov 
boundary conditions'^ for g{x) hold. 
We obtain the noise correlations from 
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Here we have J(x) = — l/(8e) J deTT[fKJ{x)] and 
Sli = li — li is the deviation of the current from its 
quantum mechanical expectation value. In this article 
we take j = 3 and thus with z = 3 we get the noise 
5*33 = S and with z = 1, 2 the cross correlations. The ef- 
fect of the current- voltage characteristics on the current 
fluctuations may be eliminated by considering the effec- 
tive charge g'eff = (3/2)|d5'/d73|, where the factor 3 arises 
from the diffusive nature of the transport. The effective 
charge yields information on the charge transferred and 
also on the energy-dependent correlations between charge 
transfers in the transport process. The matrix current in 
the first order in x 

J^^^(x) = -2id^J{x)\^^o = -LGoigodxgi ^ gidxgo) 

(3) 

is defined so that the Usadel equation in the first order 
in X is identical to Eq. ^ with the substitution g 
^1, J ^ J^^\ The Nazarov boundary conditions for J^^^ 
are given hy-^^ 

j(i) = _2Gb£^^, i = [4 + T4{^o, gs} - 

B = 4.[gi,gs] + 2Tn{gsgogigs - Mi - [gi^gs])- (4) 



Here {T^} are the eigenvalues of the transmission ma- 
trix through the interface, with conductance Gb = 
Tn/(7rfi). Below, we assume a transparent con- 
tact, Gb ^ Gd' The normalization of g{x) implies 
{^o(^), ^i(^)} = 0. This is readily satisfied by intro- 
ducing the change of the variables gi{x) = [go{x)^ ^{x)]. 

We find a parametrization for ^{x) valid also in the 
presence of a supercurrent: 
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withn = rii+ri2i, = r3i+r32i, andrn, ri2, rsi, r32, 
/co, G M. With this parametrization, x has to be gener- 
ated in the normal terminal and an arbitrary number of 
superconducting terminals be at zero potential. Because 
of the specific matrix structure of the Usadel equation, /, 
f, a, and k may be solved consecutively, and r and d are 
related by the retarded-advanced symmetry. At the NS 
interface, Eq. @ yields the boundary conditions 
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while at the normal-terminal interface the boundary con- 
ditions read gi = [fK^gN,o]^ and one may choose, e.g., 
ri=rs = ko = 0, ks = -1. 

We obtain two differential equation systems which can 
be solved consequtively, one for the retarded part (upper 
left 2x2 matrix) and one for the Keldysh part (upper 
right 2x2 matrix) of Eq. ([T]) in the first order in x- Not 
all the coefficients in the retarded part are independent 
but the equations take the form 



'^11 '11 



-7^: 



'^31 '11 

_'r)(2) // 
'^32 '11 

>(0) 



" ^^11 ^11 



^^12 ^12 ^ ''^12 ^12 



(0)^ 
-12 ^12 


I 7^(2) // 

^ ''^31 ^31 


+ '"^?>1 ^31 " 


1 77(2) // 

^ /<'32 ^32 ~ 


H7e« 


^32 


(2) // 
12 '11 


'^12 '11 


-<rn- 


I 77(2) // 








7^(2) // 
'^32 '31 


~ '^L^^si " 


I 77(2) // 




^32 



' I 11 ^11 



+Kiri2 - V- 



(2)„// 

31 '31 

7r>(l)„/ 
' 12 '11 

(2)„// 

32 '31 



-n^?Vii 

^31 ^31 
-^}^11 
^32 ^31 



77(2) // \Ty\ 
^^32 ^12 ^ ' 12 '12 



(2)„// 
32 '32 " 

'^31 ^12 



(2)„// 
31 '32 



~ '11 ^12 

' ^31^^32 = ^4- 



(8) 



Ci G M depend on 6>, ^, /l,t and their 
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Here 1Z^j , V^j ,Ui t aepena on u,(p,jL,T 
derivatives. The first and second lines of Eq. (|8|) are 
obtained by operating with Re[Tr(fi-)] and Im[Tr(fi-)] 
on the retarded part of Eq. ([1]), respectively, while 
Re[Tr(f3-)] and Im[Tr(f3-)] yield the equations on the 
third and fourth lines. The full expressions, however, 
are too long to write here. The Keldysh part obeys two 
coupled differential equations 
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Here /Cq^s^^ , Qq^s^^ ' ^1,2 ^ ^ depend on 0, 0, /l,t, ^1,3 and 
their derivatives. The first and second line of Eq. ([9|) are 
obtained from the Keldysh part of Eq. ([T]) in the first 
order in x by taking the traces Re[Tr(-)] and Re[Tr(f3-)], 
respectively. 

Putting all together, the spectral equations for 6>, (j) 
and the kinetic equations for /l^/t are first solved, 
then Eq. ([8]) for ri,r3, and thereafter Eq. ^ for 
ko^ks. Equations ([2|) and ([3|) yield a lengthy expres- 
sion for noise correlations in terms of the parameters 
^, ^, /l, /t, ^1, ^3, ^0, ^3 into which the values of these pa- 
rameters are finally substituted. Essentially because of 
the finite "coherence" parameter ri, noise deviates from 
its incoherent value^, and due to a finite supercurrent, 
rs may be finite. However, in up- down (1-2) symmetric 
structures, rs and ko always vanish in the control wire. 
We have developed a computer code to solve numerically 
these equations and present the results below. We first 
consider up-down symmetric setups and then the influ- 
ence of breaking this symmetry. 

The full phase and voltage dependence of ^eff with 
^1,2,3 = is illustrated in Fig. [2l The I — V char- 
acteristics in such a structure may be calculated from 
J{x = 0) in a way explained, e.g., in Ref. 20. Our results 
for Acj) = 0, TT coincide with those in Ref. 4 and we obtain 
a minimum of g'eff at about A(j) = O.GSn (correspond- 
ing to the maximum of the spectral supercurrent^^), 
eV = 0.5 Et. The nonmonotonic voltage dependence 
of qef[ may be understood by studying the Andreev re- 
flection eigenvalue density of a diffusive wire, which at 
V = takes the form of the Dorokhov distribution*^ 
The behavior of the noise parameters as a function of 
suggests that the returning of geff to 2e at low volt- 
ages may also be attributed to the depression of the lo- 
cal density of states at low energies. ""^^ This is illustrated 
by the fact that the voltage at which the minimum of 
g'eff is obtained follows the phase-dependent minigap in 
the superconductor-normal metal-superconductor (SNS) 
system. 

In Fig. [31 g'eff vs eV is plotted for a vanishing and 
for different cross sections Ac and lengths L3 of the con- 
trol wire, with given Et = hD/L^. The influence of a 
finite is exemplified in the inset, where the results for 
the same parameters as in the Fig.O but for A<p = 27r/3, 
have been plotted. The behavior of the supercurrent in 
somewhat similar situations has been studied earlier in 
Refs. [25I and [H. For Ac ^ 0, the voltage at which 
the minimum geff is obtained exactly coincides with the 
phase-dependent minigap in the SNS system. Generally, 
enlarging the width of the control wire or varying the 
lengths away from the symmetric case ^^1,2,3 = L tends 
to make the dip shallower. With Li^2 = L and L3 smaller 
(larger) than L, the minimum in ^eff is shifted to higher 
(lower) voltages. This is in agreement with the conclu- 
sion that the dip is caused by the anticorrelation between 
subsequent Andreev pairs as the correlations of the pairs 
manifest themselves at the length scale ~ (HD / eY^'^ . 

Under the sign reversal of A^, at energies below eV ^ 




FIG. 2: (Color online): Effective charge ^eff vs voltage and 
A(/). The voltage at which the minimum ^eff is obtained fol- 
lows the phase-dependent minigap in the SNS system (bold 
dashed black curve). Note that the junction undergoes a 
TT-transition, in which the supercurrent changes its sign, at 
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FIG. 3: (Color online): Effective charge ^eff vs eV for A(/) = 
(main figure) and A(/) = 27r/3 (inset). The cross section Ac is 
close to (blue dotted), A/4 (blue solid), A (black bold o), 
\A (red dashed) and for these curves ivi,2,3 = ^- The values 
for L3 are 1/8 (green dash-dotted □), 1/2 (green solid *), 2 
(magenta solid \j\ 4 (magenta dashed +) times Li,2 = 1^ 
with Ac = A] Et = hD/L^. For Acj) = 27r/3, the minima of 
^eff occur at lower V and except for the curves for Ac — 
and L3 = L/8, the dips are deeper than for Acj) = 0. Note 
the different energy scales on the x-axes in the inset and in 
the main figure. 



the dissipative and superconducting parts of the spectral 
charge current jt remain invariant i2i A direct calculation 
shows that also Sij is invariant under the sign reversal 
of A(j) (also in asymmetric structures). This means that 
the quasiparticle current is in no way correlated with the 
supercurrent flowing in the system, although the pres- 
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FIG. 4: (Color online): Effective charge ^eff for different val- 
ues of phase difference in left-right symmetric (s = 1) and 
asymmetric (s = 1/4) structures. 




FIG. 5: (Color online): Normalized differential cross correla- 
tions (l/2e)(d5'23/d/3) vs eV , for A(/) = (black solid), A(/) = 
7r/2 (blue dotted), A(/) = 27r/3 (red dashed), A(/) = 87r/9 
(green dash-dotted) and s — 1/4,1/2,3/4,1 for the curves 
from top to bottom. The cross correlations related to the op- 
posite wire are directly obtained from S\z — —S23 — Sss- The 
results for 5*23 and the symmetry parameter s are identical 
for those for S13 and 2 — s by the symmetry of Sij under the 
reversal of the phase gradient. 

ence of the latter changes the correlations in the pre- 
vious. Hence in the left-right symmetric structures, we 
have Sis = S23 = —S33/2 through current conservation. 
Note that if the supercurrent would be replaced by a 



dissipative current, the cross correlations would depend 
on the relative signs between Is and the "circulating" 
current Ii — h- In order to study the effect of asym- 
metry, we introduce a symmetry parameter, < 5 < 2, 
measuring the distance between the wire 3 and Ti such 
that Li/ s = L2/ {2 — s) = Ls = L. The symmetric sys- 
tem described above corresponds thus to 5 = 1. The 
same results for g'eff apply for the symmetry parameters 
s and 2 — 5 as ^eff is invariant under the change of the 
sign of A(j). Hence we restrict ourselves to s < 1 with- 
out loss of generality. In Fig. [H we have calculated g'eff, 
i.e., a normalized autocorrelation function by varying s 
and A^. With decreasing 5, the coherence in the control 
wire increases as the other superconducting terminal is 
brought closer to it. At voltages higher than Er/e^ but 
in the region where the coherence is not fully suppressed, 
the enhanced coherence gives rise to the long tails with 
(7eff < 2e in Fig. [H However, at voltages near the mini- 
mum of (7eff, decreasing s suppresses the dip. 

Figure [5] represents the normalized differential cross 
correlations (l/2e)(d5'23/d73). In the lowest curves for 
the symmetrical case, s = 1, we have (l/2e)(d5'23/d/3) = 
— geff/6e. In the incoherent region eV ^ Et^ the abso- 
lute value of 6*23 diminishes linearly with decreasing 5, as 
one may anticipate on the basis of the Kirchoff rules. In 
the coherent regime, d6'23/d/3 reflects, e.g., the relative 
changes in d/2/dF and d/3/dV with V. In the symmet- 
ric case, these relative changes have equal magnitudes. 
If the former exhibits a larger change than the latter, 
(l/2e)(d5'23/d/3) may obtain larger negative values than 
in the incoherent regime. With a finite supercurrent these 
dips correspond to the processes in which two electrons 
are injected from the normal reservoir and an Andreev 
pair enters the superconductor T2. 

In conclusion, we have found a physically transparent 
and computationally efficient way to calculate the full 
phase and voltage dependence of the noise correlations 
in mesoscopic diffusive wires in the presence of supercur- 
rent. We found that the strength of the anticorrelations 
between the Andreev pairs flowing in the structure is 
closely related to the magnitude (but not to the sign) of 
the spectral supercurrent and the variations in the local 
density of states. 
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